home
***
CD-ROM
|
disk
|
FTP
|
other
***
search
/
Turnbull China Bikeride
/
Turnbull China Bikeride - Disc 1.iso
/
ARGONET
/
PD
/
MATHS
/
RLAB
/
RLAB125.ZIP
/
!RLaB
/
rlib
/
tril
< prev
next >
Wrap
Text File
|
1994-09-23
|
812b
|
41 lines
//-------------------------------------------------------------------//
// Syntax: tril ( A )
// tril ( A , K )
// Description:
// tril(x) returns the lower triangular part of A.
// tril(A,K) returns the elements on and below the K-th diagonal of
// A.
// K = 0: main diagonal
// K > 0: above the main diag.
// K < 0: below the main diag.
// See Also: triu
//-------------------------------------------------------------------//
tril = function(x, k)
{
if (!exist (k)) { k = 0; }
nr = x.nr; nc = x.nc;
if(k > 0)
{
if (k > (nc - 1)) { error ("tril: invalid value for k"); }
else
if (abs (k) > (nr - 1)) { error ("tril: invalid value for k"); }
}
y = zeros(nr, nc);
for(i in max( [1,1-k] ):nr)
{
j = 1:min( [nc, i+k] );
y[i;j] = x[i;j];
}
return y;
};